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ABSTRACT 

Seyfert narrow line region (NLR) emission line ratios are remarkably uniform, 
displaying only ~0.5 dex variation between galaxies, and even less within an in- 
dividual object. Previous photoionization and shock models of this region were 
unable to explain this observation without the introduction of arbitrary assump- 
tions or additional parameters. Dusty, radiation pressure dominated photoion- 
ization models provide a simple physical mechanism which can reproduce this 
spectral uniformity between different objects. In the first paper of this series 
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we described this model and its implementation in detail, as well as presenting 
grids of model emission lines and examining the model structures. Here we ex- 
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plore these models further, demonstrating their ability to reproduce the observed 
Seyfert line ratios on standard line diagnostic diagrams in both the optical and 
UV. We also investigate the effects that the variation of metallicity, density and 
ionizing spectrum have upon both the new paradigm and the standard photoion- 
ization models used hitherto. Along with the standard diagnostic diagrams we 
provide several new diagnostic diagrams in the UV, Optical and IR. These new 
diagrams can provide further tests of the dusty, radiation pressure photoioniza- 
tion paradigm as well as being used as diagnostics of the metallicity, density and 
ionizing spectrum of the emission line clouds. 

Subject headings: galaxies: active — galaxies: Seyfert — ISM: general — line: 
formation 



1. Introduction 



The emission lines of active galaxies have often been used in conjunction with models 
to constrain the physical and ionization structure of the emitting regions. In particular 
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ratio diagrams or line diagnostic diagrams prove to be an excellent visual aid in interpreting 
the emission line data. First used systematically by Baldwin, Phillips, & Terlevich (1981), 
line diagnostic diagrams take the observed ratios from emission line galaxies and create a 
two-dimensional classification scheme which is better able to differentiate between excitation 
mechanisms and determine other properties of emitting regions such as density or chemical 
abundance in the gas phase. The line diagnostic diagrams by Veilleux & Osterbrock (1987) 
are particularly useful as they involve ratios of lines which are not greatly separated in 
wavelength and so minimize the effects of differential reddening by dust within the emission 
line region. These diagrams are capable of distinguishing three different groups of emission 
line galaxies: those excited by starbursts and two classes of object excited by an active 
nucleus - the Seyfert narrow line regions (NLRs) and the low ionization nuclear emission-line 
regions (LINERs). These diagrams are additionally interesting in that they show that the 
emission from NLRs is remarkably uniform, with only ~0.5 dex variation between Seyferts 
and less within individual galaxies. This uniformity of the spectral properties has since been 
confirmed in much larger samples (eg Veron-Cetty & Veron 2000). 

The use of line diagnostic diagrams (LDDs) is not limited to distinguishing the excitation 
source of active galaxies. They can also help us understand the details of the physical 
processes going on within the ionized nebulae by comparing observations with the predictions 
of both photoionization and shock models (Evans & Dopita 1985; Dopita & Sutherland 
1995). The best fitting models can, in turn, indicate the physical parameters and excitation 
mechanisms of the NLR. 

Such comparisons have led to the acceptance of both of the mechanisms for excitation 
in narrow line regions. However, excitation by fast, radiative shocks (Dopita & Sutherland 
1995, 1996) appears mostly in the extended NLR (ENLR) associated with radio galaxies and 
in LINERs, and is applicable only for a few Seyfert galaxies. The photoionization models, 
and variations thereof (e. g. Komossa & Schulz 1997; Binette et al. 1996, 1997; Baldwin, 
Ferland, Korista, & Verner 1995; Ferguson et al. 1997), are able to reproduce the Seyfert 
observations with only a few failings. The main problem with the standard photoionization 
models is that they are unable to provide the observed uniformity in emission line ratios 
without making arbitrary (and possibly unphysical) assumptions. The less than 0.5 dex 
variation in observed line ratios requires an approximately constant ionization parameter of 
U ~ 10 -2 , where the ionization parameter is a measure of the number of ionizing photons 
against the hydrogen density (U = S^/nnc). Since there is no physical reason why the ratio 
of photon density to gas density should be constant within a given Seyfert, or from one 
Seyfert to another, such a result is puzzling to say the least. Furthermore, many Seyferts 
show strong coronal lines of highly ionized species which can only be produced in a plasma 
with U ^1.0. So why do low and intermediate excitation species indicate U ~ 10~ 2 ? 
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In order to account for these failings of the standard models we have proposed a new 
paradigm for the photoionization of the NLR clouds, that of dusty, radiation pressure dom- 
inated photoionization (Dopita et al. 2002, hereinafter DG02). 

We have demonstrated previously that this model provides a self-consistent explana- 
tion for the remarkable similarity between the emission line spectra of NLR. The radiation 
pressure acting upon the dust provides a simple controlling factor for the moderation of the 
density, excitation and surface brightness of the emission line region. This limits the models 
of the low- and intermediate-excitation line ratios at high ionization parameters to ~0.5 dex 
in the line diagnostic diagrams, but allows the coexistence of high ionization parameter and 
low density regions within a single model. 

Following the introduction of this model, we gave details on how this model was imple- 
mented within the photoionization & shock code MAPPINGS ill and a grid of UV, Optical 
and IR emission lines, covering a range in metallicity, density, power-law ionizing continuum 
and photoionization parameter (Groves, Dopita & Sutherland 2003, Paper 1 of this series). 
We also examined the physical structure of the dusty model and how this varied with the 
different input parameters. 

We continue this work here, presenting line diagnostic diagrams for both our dusty 
model and the standard photoionization model derived from the line ratios given in Paper 
1. 



2. Dusty, Radiation Pressure Dominated Models 

The inclusion of dust into photoionization models affects the final emission spectrum in 
several ways. As well as simply absorbing EUV radiation and competing with Hydrogen for 
the ionizing photons, dust affects the temperature structure of the NLR clouds through the 
process of photoelectric heating. 

In order to be physical, an isobaric photoionization model must include the effects of 
radiation pressure. The force of radiation can be imparted to both the gaseous medium 
and dust, and results in a radiation pressure gradient. Since the grains are charged and are 
therefore locked to the plasma by coulomb forces, the radiation pressure gradient on dust 
results in a gas pressure gradient identical in size. Standard photoionization models are 
isochoric and therefore cannot take this effect into account. To demonstrate these effects 
we have run the standard isochoric model and the new dusty, radiation pressure dominated 
model over a large set of input parameters. 
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With simple calculations it is easy to show that at an ionization parameter of log U ~ 
—2, dust begins to dominate the opacity of the ionized cloud and hence the radiation pressure. 
It is also around this value of the ionization parameter that radiation pressure developed at 
the ionization front of the NLR cloud becomes comparable with the gas pressure. Therefore 
at high ionization parameters, the pressure in the ionized gas, and hence density (since the 
electron temperature is always ~ 10 4 K), is determined by the external ionization parameter, 
Uq and the local ionization parameter becomes independent of the external ionizing flux. 
The result of this is that at high ionization parameter the emission line spectrum of the 
low- and intermediate ionization species is effectively independent of the external ionization 
parameter. This independence was illustrated through the example line diagnostic diagrams 
(LDDs) given in DG02 and is demonstrated again in this paper. 

Once the exciting mechanism of the NLR of the active galaxies has been clearly iden- 
tified, the diagnostic line ratio plots will also be useful in defining the physical conditions 
within the NLR. To facilitate this, we have run a series of the dusty, radiation pressure 
dominated models covering a range in density, metallicity, power-law index of the ionizing 
spectrum, and photoionization parameter. To provide a comparison we also present a series 
of dust-free, standard photoionization models covering the same range of parameters. We 
described the implementation of these models in depth in Paper 1 of this series (Groves, 
Dopita & Sutherland 2003), where we gave the emission line strengths from these models in 
a sequence of tables covering the parameter space. These tables were used to construct the 
diagnostic plots presented here. 



2.1. Parameter Space 

The parameters were chosen to cover the range of values that would be reasonably 
expected to be found in the NLR of Seyfert galaxies. The spread in parameter space allowed 
an examination of how each parameter affected the models, and was great enough that the 
resulting line emission from each model was distinct. 

Three Hydrogen number densities (iih) were modeled; nH = 10 2 , 10 3 and 10 4 cm~ 3 . For 
the standard isochoric (constant density) models such modeling is straightforward. However 
for the isobaric, dusty, radiation pressure dominated models the concept of a single density 
for the nebula no longer holds. Thus we have used the region near the ionization front, where 
nnn/nni ~ 1.0, as the fiducial point at which to set the density. It is in this region that the 
density sensitive lines like [S ll]AA6717, 30A or [O ll]AA3727, 9A are at their strongest and so 
this region defines the effective electron density that line observations will measure. 
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Five Metallicities were examined in the models; O.25Z , O.5Z , 1Z , 2Z and 4Z . The 
abundance set adopted for solar metallicity is given in table 1. The individual abundances of 
most elements scales with the metallicity, X/H = Zx (X/H) Q . The two contrary elements 
are Helium and Nitrogen. For Helium, the chemical yield from stars adds only a small 
amount to the primordial abundance, 

He/H = 0.0737 + O.O293Z/Z (1) 

The nucleosynthetic status of Nitrogen is unusual in that it has both primary and 
secondary nucleosynthetic components. The Nitrogen abundance at different metallicities is 
given by; 

[N/H] = [O/H] (10- L6 + i ( 2 - 33 + lo Sio[°/H])) . ( 2 ) 

The origin of this and equation 1, as well as the justification for the Solar abundance set are 
discussed in Paper I. 

Because dust is present, the gas-phase abundances are depleted in comparison to the 
total abundances. The depletion fractions adopted are given in Table 2. In the new radiation 
pressure dominated models the metals are depleted onto the dust, but for the dust-free 
standard models this depletion is artificial, with the metals effectively lost. This is done 
to ensure that the gas phase abundances in both sets of models are the same, because the 
heavy elements in the gas phase determine the cooling to a large extent and also have some 
effect upon the ionization state within the nebula. 

We use a simple power-law to represent the spectral energy distribution (SED) of the 
ionizing source, with 

F v (X V a IVm < V < ^max- (3) 

and z/ m i n = 5eV and z/ max = lOOOeV. We investigated four values of the power-law index a, 
-1.2, -1.4, -1.7 and -2.0. These indices encompass the 'standard' values usually adopted for 
modeling of the AGN spectrum. The factor of proportionality, which determines the total 
radiative flux entering the photoionized cloud is set by the ionization parameter at the front 
of the cloud, U = S+ffooc), where S* is the entering flux of ionizing photons, n the initial 
density and c the speed of light. 

This ionization parameter is the final parameter allowed to vary in this family of models, 
with log£/o varying between 0.0 and -4.0, sampled at intervals of -0.3, -0.6 and -1.0 dex. 
For the standard isochoric models the initial density is a known quantity and hence the 
ionization parameter at the front of the cloud is easy to determine. In the dusty, isobaric 
models the initial density is not known a priori. To estimate this value we assume a front 
end temperature of T = 20, 000 K and obtain the density from the set P . Depending 
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upon the temperature reached at the front of the model, the true density will vary, being 
overestimated at higher temperatures and underestimated for cooler temperatures. This 
means that the true ionization parameters in the upper range (logC/o > —2.0) are actually 
larger than that given in the text, and are smaller for the lower range (logC/o < —2.0). The 
general shape of the curves in the line diagnostic diagrams is still correct. However, the true 
range of ionization parameters is somewhat larger than indicated on these diagrams. 

3. Diagnostic Diagrams 

There are several considerations we must take into account in the choice of line ratios 
to use in the line diagnostic diagrams (LDDs). These are: 

1. The lines must be strong enough to be measurable in most Seyferts. 

2. The line diagnostic diagram must provide physical insight into the nature of the emit- 
ting system. For example, the LDDs put forward by Veilleux & Osterbrock (1987) 
provided clear diagnostics of the excitation mechanism, with NLRs, LINERs and H II 
galaxies well separated. Another good example are the UV diagnostics of Allen, Dopita, 
& Tsvetanov (1998) which are able to discriminate between photoionization models, 
shock models, and shock+precursor models. 

3. To minimize wavelength sensitive effects such as reddening and flux calibration, the 
wavelength separation between the lines that make up the ratio should be as small as 
possible. For reddening this becomes more important at shorter wavelengths, especially 
around the 2200A extinction feature seen in our galaxy. 

4. The lines should be all observable using a common instrument and the same technology. 
For example, a UV line ratio versus a Visible line ratio may be a good diagnostic, but it 
is highly unlikely that one could ensure the same aperture for these observations. Hence 
line ratio diagrams should be made up of ratios that lie within a single wavelength 
regime easily accessible by the same instrument. 

5. Emission lines which are usually blended with other lines must be either avoided or 
taken as a blend, as deblending lines can greatly increase the uncertainty of the mea- 
sured flux. 

6. One should avoid the use of lines which carry large flux errors resulting from the theo- 
retical models used. An example of this is [O l]A6300A which, though easily measurable, 
is highly dependent upon the point at which the model has been truncated. 
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3.1. Standard Optical Diagnostic Diagrams 

The line diagnostic diagrams put forward by Veilleux & Osterbrock (1987) (V&O) have 
become standard diagnostics for emission line galaxies. This is because the diagrams not 
only consist of the some of the strongest emission lines which are easily accessible to ground 
base telescopes, but they also clearly distinguish between the excitation mechanisms of the 
emission line galaxies. 

Figure 1 shows the first of the V&O diagrams; [N ll]A6583/Ha versus [O IIIJA5007/H/3. 
Figure la demonstrates the effect of different metallicities upon the dusty models, with the 
Seyfert 2 observations from the original V&O paper given for comparison. The models 
plotted are for a density of 1000 cm -3 and a power law index of a = —1.4, and cover a range 
in ionization parameter of —4.0 < logf/ < 0.0. Lines of constant ionization parameter are 
marked and labelled at set intervals. Each model curve represents a different metallicity, 
increasing from left to right, with the metallicity of each curve labelled. Also indicated on 
the diagram is the effect of reddening upon the models. The direction and magnitude of the 
arrow indicates the effect of an external dust screen with 10 magnitudes of visual extinction 
using the reddening curve from Calzetti (2001) (their equation 8). The diagram reveals the 
large spread in the [Nil] /Ha ratio for the set of photoionization curves, which is largely due 
to the varying contributions of the secondary component of Nitrogen nucleosynthesis. A 
comparison of the photoionization curves with the observational data also reveals that most 
of the NLR in Seyfert galaxies are characterized by super-solar metallicity, with a mean value 
of about 2 Z Q . The sub-solar curves are unable to reproduce the data. 

Using this information, figure lb compares both the dusty models and the standard, 
dust-free isochoric models with the observational data. The different sets of models are 
labelled and cover a range of power-law index from a = —1.2 to —2.0. Both sets of mod- 
els have a metallicity of 2 Z Q and nH=1000 cm -3 , with each model curve ranging from 
logf/o = —3.0 to 0.0. For both sets of models, as the power-law becomes flatter (i.e. harder 
spectral energy distribution), the trajectories of the models are displaced from the bottom- 
left towards the top-right region of the diagram. The dusty models "stagnate" in the region 
occupied by the observational data at high ionization parameter. That is, as the ionization 
parameter increases, the dusty models tend to return the same values for the ratios, and be- 
come degenerate in terms of U at the highest values. As for the dust-free models, [Nil] /Ha 
decreases continually with increasing ionization parameter, as the overall ionization state of 
the gas increases. This behavior provides a clear distinction between the models and is one 
of the principal reasons successes of the dusty, radiation pressure dominated photoionization 
models, as discussed in Dopita et al. (2002). 

There are several other points to note within the diagram. One is the difference between 
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the behavior of the two sets of models in the [O lll]/H/3 ratio at low ionization parameter. At 
high ionization parameter a difference is expected due to the effects of extinction by dust and 
radiation pressure on dust, but at low ionization parameter such effects should be minimal 
as dust is no longer the dominant opacity. The difference between the models arises due to 
the definition of the ionization parameter at the front of the cloud, as discussed in §2.1. In 
the case of log£/o of -3.0, the dusty models overestimate the effective ionization parameter 
due to the overestimate of the assumed initial temperature. Hence the dusty models display 
a weaker [O ill] /Ha ratio than the dust-free set. The apparent difference is therefore an 
artifact of the initial condition, and should be ignored. The lack of observational points 
below log U ~ —3 is not a failing of the models, rather it is due to the definition of Seyfert 2 
galaxies. The region below that occupied by the Seyferts is not empty, rather it contains the 
class of galaxies classified as LINERs (see eg. V&O). LINERs may include a subset which 
could be considered low power Seyferts, but they also include shock excited emission line 
galaxies and possible old starbursts as well, and thus we have not considered these within 
our sample. 

The other point to note is the spread in model sets due to the variation in power-law 
index, a. For both sets of models a flatter (harder) ionizing spectrum leads to an increase 
in both ratios, especially [Olll]/H/3. This is due to the increase in photons available to 
ionize these higher ionization species compared to hydrogen. At low ionization parameter, 
the region covered on the diagram by the models of different power-law index is similar to 
that found for models with different metallicity. This indicates the limited utility of these 
ratios as either a Z or an a diagnostic. At high ionization parameter however, the variation 
of a is clearly distinguished from that of Z in the dusty models, with the models at different 
a curving in upon each other such that they deliver little spread in the [Oin]/H/3 ratio or 
the [Nil] /Ha ratio. The observational data is well covered by the models, and the most 
reasonable values of a probably lie in the range —1.4 < a < —1.2. 

Figure lc shows the V&O diagram [O l]A6300/Ha versus [O lll]A5007/H/3, with the data 
points being taken again from the same paper. Though the strength of the [Ol]A6300 line 
may not be entirely trustworthy for the reasons discussed previously, the line is strong and 
easily measured and the ratio proves to be a good diagnostic. As discussed in Paper 1, 
our models are truncated at the point at which Hn/H ratio drops below 1%, such that 
the final temperature is a few 100 K. This means that the majority of the O I emitting 
region is encompassed within the model as the emission of O I drops below 1000K. Figure lc 
demonstrates the effect of variation of metallicity on both the dusty and the dust-free models. 
The models have Uq ranging from 10~ 3 to 10°, with a fixed at —1.4 and the density at 1000 
cm~ 3 . The metallicity increases from 0.25Z Q for the leftmost curve to 4Z Q on the right. As 
both line ratios are ratios of oxygen to hydrogen lines, any variations are due to changes in 
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the temperature and ionization of the gas. The multi- valued nature of the curves arises due 
to the combination of two effects. As metallicity increases, the strength of the metal emission 
lines initially increases due to the increase in abundance. However, increasing the metallicity 
also increases the cooling efficiency of the gas and hence lowering the temperature. At high 
ionization parameter, there is a clear distinction between the dusty and dust-free models, 
as seen in the previous diagrams. The dusty model line ratios stagnate in a limited zone of 
parameter space, while the dust free models display a large decrease in the [O i]/Ha ratio with 
increasing ionization parameter. As in the previous diagnostics, the observed data points are 
best fit by the dusty models, with the curves at high ionization all clustered in the region of 
interest. The data is best reproduced by the curves with a metallicity of Z — 1 — 2Z &) but, 
because oxygen is a primary rather than secondary nucleosynthesis element, the distinction 
between the curves of different metallicity is not as clear as in the case of the [Nil] /Ha ratio. 

Figure Id shows the effects of variation of a upon the dusty model ratios. The models 
all have Z = 2Z Q , % = 1000 cm~ 3 , and each curve covers the range —3.0 < log[/o < 0.0. 
Each curve represents a different power-law index, as labelled. This diagram more clearly 
demonstrates the spread due to a at low metallicities, and the convergence of the models in 
a narrow strip of parameter space, similar to that seen in figure lb. Though not shown here, 
the dust-free models also appear very similar to the dust-free curves of figure lb as well. 
This diagram clearly indicates how well the dusty models reproduce both the clustering and 
the absolute values of the observational data. A value of a between —1.4 and —1.9 provides 
the best fit to the observed range of observational points. 

The final V&O diagram is [S ll]AA6717, 30/Ha versus [O m]A5007/H/3. Figure le demon- 
strates the effects of metallicity variation upon the two sets of models, while figure If demon- 
strates the effects of power-law index variation upon both. Both figures use the same set of 
fiducial parameters as before. 

In most ways these diagrams are similar to the previous figures. They display a large 
spread due to metallicity and power-law index, and show similar shapes to the previous 
model curves. At high ionization parameter, the dusty and dust-free models separate as in 
the previous figures, with the dusty models displaying the same stagnation in the region 
of interest. The observational data is similarly best reproduced by the dusty models, with 

a metallicity Z ~ 1Z & and a power-law index a ~ —1.4 1.7. Note that though the 

Z ~ 1Z & seems to better fit the data when metallicity variation is considered, the diagram 
of the a variation (using Z = 2Z Q to remain consistent) demonstrates that the variation in 
power-law index may also be quite large. 

Of the three V&O diagnostic diagrams, the [N n]A5683/Ho; is the clearest metallicity 
diagnostic thanks to the effects of the secondary nucleosynthetic component of the nitrogen 
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production. It has well separated metallicity curves that cover a range which is much broader 
than that due to power-law index variations. This means that it is not degenerate with a 
as the other two diagrams are. This diagnostic diagram shows that the mean metallicity 
in Seyfert NLRs is ~ 2Z Q although it may truly range from 1Z Q up to 4Z . The result of 
a higher metallicity is not surprising, as the new abundance set used and the depletion by 
dust both act to reduce the gas abundance compared to previous models. A value of ~ 2Z Q 
is approximately the same as the old Solar values (Anders & Grevesse 1989), as has been 
estimated previously for NLR metallicity. 

An estimate of the average power-law index can be obtained by consideration of all 
three diagrams. In these, the range —2.0 < a < —1.2 encompasses all of the observed data, 
but more typical values of the index range between —1.4 and —1.7. 

The mean density of the Seyfert NLR clouds is harder to determine as all three diagnostic 
diagrams show small variation due to density at high ionization parameters when compared 
to the effects of the other parameters, especially metallicity. 

The self-consistency of the modelling is demonstrated by the fact that the observed ratios 
on all three diagnostic diagrams are reasonably reproduced by the same set of parameters, 
with a very broad range of ionization parameters able to reproduce the tight clustering of 
the data. 



3.2. Standard UV Diagnostic Diagrams 

Whilst the optical diagnostic diagrams provide an excellent tool for separating star- 
burst galaxies from those with an AGN, these optical line ratios are unable to conclusively 
distinguish between shocks and photoionization as an excitation mechanism for active galax- 
ies. Ultraviolet emission lines, however, tend to be much stronger in shocks than in simple 
photoionization models and hence much better diagnostics for shock and photoionization sep- 
aration. This difference is due to the collisionally excited lines in the UV, like Civ A 1549, 
having increased emission at the higher temperatures found in post-shock gas (Dopita & 
Sutherland 1995, 1996, and references therein). This temperature sensitivity also proves to 
be useful in distinguishing between classical photoionization models and the dusty models 
presented in this paper. 

A further benefit of UV diagnostic diagrams is that in high redshift active galaxies, the 
UV emission lines are shifted to the optical band, and can be observed with ground-based 
telescopes. Hence the UV diagnostic diagrams are better for determining the excitation 
mechanism of the high-z galaxies than the standard optical ones which may be difficult or 
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impossible to observe. 

UV diagnostic diagrams have been developed and explored in the case of high-Z radio 
galaxies by Villar-Martin, Tadhunter, & Clark (1997) and comparisons of shock and pho- 
toionization models have been made by Allen, Dopita, & Tsvetanov (1998) (hereafter ADT). 
Here, we continue this exploration and comparison with the inclusion of dusty photoioniza- 
tion models. 

The first three UV diagrams of Villar-Martin, Tadhunter, & Clark (1997) use three 
UV line ratios; Civ A1549/C lll]A1909, Civ A1549/He IIA1640 and Clll]A1909/HellA1640. 
The observational data on each figure comes from ADT, with the stars marking the high-z 
radio galaxy (HZRG) observations and the triangles representing HST observations of three 
nearby Seyferts NGC 1068, NGC 5643 and NGC 5728. In addition to the observational data 
each figure is marked by an arrow indicating the size and direction of the reddening vector 
by a foreground dusty screen with an Ay of three magnitudes. As for the previous section, 
the Calzetti (2001) curve is used for reddening correction. 

Figure 2 shows the effect of metallicity variation on the Civ A1549/Clll]A1909 versus 
Civ A1549/Hell A1640 diagnostic diagram. The models all have a power-law index of 
a = —1.4, and density nn = 1000 cm -3 . Each curve covers a range in ionization parameter 
of —2.3 < log [To < 0.0, with each having a different metallicity as labelled. The range 
of ionization parameter shown here is smaller than in the optical diagrams as Civ rapidly 
weakens for U < 10~ 2 . As in the optical diagrams, the line ratios stagnate at high ionization 
parameter in the region where both the HZRGs and the Seyfert galaxies lie. They also 
stagnate in metallicity as both line ratios reach their maximum value at Z ~ 1Z & . Any 
value between 0.25.Z© < Z < 2.0Z Q is equally effective in reproducing the observations. A 
metallicity less than solar may even be plausible for the HZRGs which are observed at early 
stages of galaxy formation. 

Following the format of the optical diagnostic diagrams, figure 2b demonstrates the 
effects of power-law index variation upon both the dusty and dust-free models at a metallicity 
of 2Z & , and uh = 1000 cm" 3 . The index becomes steeper from top to bottom on the diagram 
for both sets of models and the ionization parameter covers the same range as before. Little 
can be inferred from this diagram other than that both the dusty and dust-free models 
provide an equally good fit to the observational points. However, it is significant that there 
are no observational points above log(C iv/ C ill]) > 0.6, just above the point where the dusty 
models stagnate. 

As discussed in Paper 1, care must be taken with these UV diagnostics due to the 
presence of resonance lines like Civ. The radiative transfer of these lines is not treated 
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exactly within the photoionization code, but rather approximated with an escape probability, 
which includes the opacity effect of dust. This means that there is a larger error associated 
with the ratios which involve these lines and they should be treated as less certain. 

The second UV diagnostic diagram is Civ A1549/C lll]A1909 versus C ill] A1909/HellA1640, 
shown in figure 2c, where the effect of metallicity variation upon these line ratios in the dusty 
models is demonstrated. Again all curves have a = —1.4 and tin = 1000 cm -3 with logC/o 
ranging from —3.0 to 0.0. The shape and form of these curves are very similar to that found 
in figure 2a. Here the effect of metallicity variation is once again more apparent at low 
ionization parameters, but both metallicity and ionization parameter become degenerate at 
high ionization parameter. All that can be said is that the dusty models fit the data points 
in the range logf/ > —2.5 and Z ^ 2.0Z Q . 

In figure 2d we assume a 2Z Q metallicity to compare the dusty models with the dust-free 
models for these ratios. Shown in this comparison is the effect of density variation on these 
two models. The power-law index is —1.4 and ionization parameter is —3.0 < \ogU < 0.0 
for these sets of models. In both sets of models the density increases from 10 2 cm~ 3 at 
the bottom to 10 4 cm~ 3 at the top. Clearly, although some density sensitivity is present, 
the spread due to metallicity or spectral index is just as large and hence these parameters 
cannot be separated by the use of these ratios alone. Nonetheless, it is clear that the dusty 
photoionization models provide a much better description of the range of observations than 
do the dust-free models. 

The third Villar-Martin, Tadhunter, & Clark (1997) diagram has C m]A1909/He IIA1640 
versus Civ A1549/He IIA1640. As before, we demonstrate first the effect of metallicity upon 
these ratios for the dusty model, in figure 2e, using a = —1.4, uh = 10 3 cm~ 3 and —3.0 < 
logf/o < 0.0. Though rotated and flipped relative to the previous diagrams, the curves on 
this diagram show a similar form to that seen on both the other two diagrams. 

Finally, figure 2f shows the effect of a variation with these ratios on both sets of models, 
using the standard Z = 2Z Q , uh = 1000 cm~ 3 and —3.0 < log£/o < 0.0. The power-law 
index increases from left to right in both sets of models. This also has the same general form 
as the previous diagrams. Once again, the line ratios of the dusty models stagnate in U and 
Z in the region occupied by the observational points. These dusty models provide a much 
better fit to the observations than the dust-free models. 

Allen, Dopita, & Tsvetanov (1998) also provided several other diagnostic diagrams, 
mainly with the purpose of distinguishing shock excitation from photoionization. These 
diagnostics also prove to be valuable in further establishing the validity of the dusty, radiation 
pressure dominated photoionization models. 
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Figure 3a shows the sensitivity to metallicity and ionization parameter of C IV A1549/C ill] A 1909 
versus Cn]A2326/Cin]A1909. The empty triangles represent the Seyfert galaxies shown in 
the previous UV diagnostic diagrams. This figure is somewhat more useful in distinguishing 
between different metallicities, but unlike the previous three UV diagnostics, it appears to 
indicate a metallicity greater than 2Z & . Another possibility is that these galaxies may also 
have some shock component to them, which, as seen in figure 2d of ADT, would bring the 
models closer to the observations. 

This diagram is also somewhat sensitive to the power-law index, as shown in figure 3b, 
which displays both the dusty and dust-free models are displayed with 2Z Q , nn = 10 3 cm -3 
and —3.0 < logC/ < 0.0. Clearly, the dusty models once again provide a better description 
of the observations than the dust-free models, which have too weak Cii]/Cm] ratio and too 
strong Civ/Cn] ratio at high ionization parameter. 

In figure 4a we examine an optical-near UV diagnostic, [Ne v]A3426/[Nelll]A3869 versus 
[Olll]A5007/H/3. The effect of metallicity variation upon these ratios in the dusty models is 
shown in figure 4a. Though the spread is small compared to some of previous UV diagnostics, 
this diagnostic has the additional benefit of the minimal error and reddening in the [O lll]/H/3 
ratio. Included on the diagram are three observational data sets. The empty triangles 
represent the Seyfert galaxies from the previous UV diagrams. The crosses are data from 
four Seyfert 2 galaxies from Allen (1998). The asterisks represent observations of Seyfert 
2 galaxies from the data set of Koski (1978). All three sets of data are consistent with 
a metallicity somewhere between 1Z Q < Z < 4Z Q . Although the dusty models stagnate 
in terms of the ionization parameter in the region occupied by the data, not as much can 
be inferred from this because the dust-free models are multi-valued in Uq in this region as 
well. This is shown in figure 4b, which plots curves of different a for the two models at 
Z = 2Z Q , n H = 1000cuT 3 and —2.6 < \ogU < 0.0. a increases with increasing [Om]/H/3. 
The two sets of models differ in their behavior in a, with the dusty model curves stagnating 
in a very restricted region of parameter space but the dust-free models becoming more 
widely separated. The fact that the stagnation point of the dusty models agrees with the 
observational points gives further credence to the general validity of this class of models. 

In addition to the ADT paper there have been more recent investigations of UV diag- 
nostics, such as Inskip et al. (2002a,b) who used the diagnostics to examine high redshift 
(z ~ 1) radio galaxies. The redshift of these galaxies brought several of the UV diagnostic 
lines into the optical range, making them accessible to ground based telescopes, with long 
enough integration times. Due to the fact that these galaxies have very energetic expanding 
radio lobes, they are expected to be excited by both shocks and by nuclear photoionization. 
Being able to distinguish between these differing excitation mechanisms is important in our 



-14- 



understanding of these high-z AGN. 

The UV diagnostic diagram chosen by Best, Rottgering, & Longair (2000b), and ex- 
panded by Inskip et al. (2002a), combines the same neon line ratio of the previous diagram; 
[Ne m] A3869/[Ne v] A3426 but plotted against the UV carbon line ratio C m]A1909/C n] A2326. 
This diagram is very good for distinguishing between shock excitation and photoionization. 
It also turns out to be very good for distinguishing between the dusty models and the dust- 
free models (figure 5a). The HZRG observations from Inskip et al. (2002a) (their figure 14) 
are marked on the diagram. Triangles indicate the 6C radio galaxies and stars the 3C radio 
galaxies. The ionization parameter is restricted to the range —2.3 < logC/o < 0.0 because 
[Ne v] is a very high ionization species and becomes very weak at lower ionization parameters. 
Figure 5a shows the sensitivity of the emission line ratios to the metallicity. Again, many of 
these radio galaxies sit in the region of the diagram that is consistent with high ionization 
parameter; log [To ~ —1.0. However, little can be inferred about the metallicity. The radio 
galaxies with weak Clll]/Cll] ratios (which are believed to be shock or shock+precursor 
excited (Best, Rottgering, & Longair 2000a,b; Inskip et al. 2002a,b)) cannot be reproduced 
by these dusty models without resorting to very high [Z > AZ Q ) metallicities. This is similar 
to figure 3. 

Figure 5b demonstrates the use of this diagram in distinguishing between dusty and 
dust-free photoionization models. The model curves, with Z = 2Z Q , nn = 1000 cm -3 , vary 
in power-law index, ranging from -1.2 at low Ciii]/Cii] to -2.0 for the models with high 
Cm]/Cii]. The dusty models are clearly able to reproduce the HZRG observations that 
are believed to be photoionized, whereas the dust-free models would imply a much flatter 
power-law (a < —1.2) or else a higher metallicity to even approach the observational points. 

Though our ability to derive a clear definitive set of parameters for the Seyfert Galaxies 
and the high redshift radio galaxies is limited in these UV diagnostic diagrams due to the 
degenerate nature of many of the curves, one detail stands out in all diagnostics; the dusty, 
radiation pressure dominated models provide undeniably the better fit to the observations 
over the standard dust-free models. In all cases the dusty models not only reproduce the 
data well, but also tend to become degenerate in terms of the ionization parameter precisely 
in the region occupied by the observations. 

3.3. Further Useful Diagnostic Diagrams 

The previous sections demonstrated the utility of the dusty, radiation pressure domi- 
nated models. We will now explore several new or relatively unexplored line ratio diagrams, 
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not only to gain a deeper understanding of the dusty photoionization models but to also 
find diagrams which can provide diagnostics of the metallicity, density and the slope of the 
ionizing power-law. 

These diagrams extend the wavelength base from the far-UV into the near and mid-IR. 
They adhere to the guidelines of Veilleux & Osterbrock (1987) re-iterated at the start of 
section 3. The program used to create these diagnostics is available at 

http://www.mso.anu.edu.au/~bgroves/linedata. Also available there is the model data 
used to create these diagnostics, which were discussed in depth in Paper 1. This includes the 
standard, dust-free, undepleted models discussed in Paper 1, if the readers wish to examine 
the diagnostic diagrams of this model. 

For the dusty models we adopt a fiducial model which has a metallicity of 2Z & , a power- 
law index of a = —1.4 and a density of = 1000 cm -3 following the results of the previous 
sections. Each of these parameters is adjusted separately to determine the sensitivity of each 
line diagnostic to each parameter in turn. For the ionization parameter, the greatest possible 
range is explored. This turns out to be quite restricted for the line ratios which involve high 
ionization species. 

3.3.1. UV Diagnostics 

The benefits of UV diagnostics in providing insight into the physics of high redshift 
active galaxies and in distinguishing between photoionization and shock excitation has been 
examined in depth in earlier works by Villar-Martin, Tadhunter, & Clark (1997); Allen, Do- 
pita, & Tsvetanov (1998); Best, Rottgering, & Longair (2000a,b) and Inskip et al. (2002a,b). 
However, with the deeper surveys of high redshift galaxies currently underway, we can expect 
that diagnostic plots involving lines in the far-UV will soon become even more useful. 

The first far-UV diagnostic takes two line ratios from ADT; Cm A977/C m]A1909 and 
NlllA991/Nlll]A1750. Both ratios are temperature sensitive and hence are good at sepa- 
rating shock excitation from photoionization (see figure 3 of ADT). When plotted against a 
ratio which increases steadily with Uq, such as Civ A1549/C ill] A1909 they also provide rea- 
sonably good metallicity diagnostics, although the spread due to a or nn are similar in range. 
Both of these line ratios have a similar dependence on the electron temperature. Thus, when 
plotted together, variations in the metallicity, power-law index and ionization parameter are 
relatively indistinguishable. However, the ratios do provide a diagnostic which is able to 
distinguish variations in the density, as shown in figure 6. This arises mainly through the 
nitrogen ratio. The N ill] A 1750 line can be represented as a five level system, containing 
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the lines that make up the 1750 line and a fine structure line at 57 jum. As the density in- 
creases, the timescale for collision from the fine structure level becomes less than the emission 
timescale, and the flux that would normally be emitted through 57 jum, is emitted at 1750 
A. As the C ill A977 resonance line is reasonably unaffected by density, this ratio becomes 
density sensitive. The density sensitivity is reversed at low U due to the lower temperature 
and hence lower collisional excitation rate. 

One problem with these ratios though is that the reddening is quite large due to the 
short wavelengths of the lines and the large wavelength separation between numerator and 
denominator. Corrections due to reddening are also quite uncertain at these short wave- 
lengths, hence the absence of the reddening arrow on figure 6. The small (~ 0.1 dex) 
separation between the three densities and the likely large absolute flux errors resulting from 
the uncertain reddening corrections means that this diagram is probably not useful as a 
density diagnostic. 

In figure 7 we plot Civ A1549/C m]A1909 versus CTllA977/NlllA991. These line ratios 
are much less sensitive to both temperature and reddening than those used in the previous 
diagram. The sensitivity to a is also minimized as a consequence of the insensitivity to 
temperature. The ratios also show very little variation due to density changes. However the 
Clll/Nlll ratio is quite sensitive to changes in metallicity. This is thanks to the secondary 
nucleosynthetic component of the nitrogen abundance. The Civ/Cm] ratio is primarily 
useful in distinguishing between different values of Uq. In addition to the effects of metallicity 
variations, what is also visible is the stagnation at high values of Uq due to dust and the 
radiation pressure upon it. 

Figure 8 exploits the metallicity dependence further by plotting Cm]A1909/Niii]A1750 
versus C IIIA977/N IIIA991. Whereas the previous diagram only separated metallicity on 
one axis, here both axes are ratios of C/N and show large separation between metallicity 
curves. In addition, the reddening corrections are small, since both pairs of lines are closely 
separated in wavelength. Thus we have a reasonably sensitive metallicity diagnostic. The 
only failing of this diagnostic is the short wavelengths of the C m/N ill ratio which may make 
this diagnostic difficult to observe. 

Our final far-UV diagnostic diagram also provides a metallicity diagnostic, plotting 
OviAA1032,8/CivA1549 versus N VA1240/C IVA1549 (figure 9). This diagnostic diagram 
should prove to be useful in several ways. Firstly, it involves only strong lines and is thus 
easily measurable. Secondly, the response of the curves to a and % variations is negligible 
compared to the metallicity variations. Thirdly, the sensitivity of the O vi/Civ to ionization 
parameter is strong. Thus this diagram provides an excellent metallicity and ionization 
parameter diagnostic for photoionized narrow line clouds. 
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3.3.2. Optical Diagnostic Diagrams 

The optical spectral range is easy to observe and provides a large number of lines from 
which to form diagnostics. Many of these line ratio diagnostics are either temperature or 
density sensitive as well. 

Figure 10 plots the near-UV ratio [Ne v]A3426/[Nelll]A3869 against the density sensi- 
tive ratio [O ll]AA3727, 9/[0 ill] A5007. We show both the density and ionization parameter 
dependence of both the dusty and dust-free models. The density increases with decreasing 
[O n]/[0 ill] ratio, from 10 2 cm -3 at the top to 10 4 cm" 3 at the bottom. Included on this di- 
agram are the observations of Seyfert 2 NLRs from Koski (1978) (asterisks) and Allen (1998) 
(crosses). Note that some of the Koski sample has been removed due to contamination by 
starbursts (see, eg Kewley et al. 2001) Within this diagram the oxygen ratio provides sensi- 
tivity to density variations, while the neon ratio provides sensitivity to ionization parameter. 
Though density variations generate a strong response in the [O ll]/[0 ill] ratio, changes in 
the other parameters also have some consequences. Variation in the power-law index affects 
the [Nev]/[Nell] ratio (as seen in figure 5a) and may be confused with the effects of density 
and ionization parameter. In terms of metallicity changes, the only major effect is seen at 
high Uq in the dusty models, where it creates a similar spread in the [O ll]/[0 ill] ratio to 
density. So, overall, this diagram provides a reasonable density diagnostic. What is obvious 
from figure 10 is that this diagram is able to strongly distinguish between the two sets of pho- 
toionization models. The large separation at high ionization parameter clearly demonstrates 
the success of the dusty model in reproducing the observations. This diagnostic ability is 
largely due to the [O n]/[0 ill] ratio and arises because of the differences in density structure 
of the isochoric, dust-free model and the isobaric, dusty model. 

The use of this ratio as a model diagnostic is explored further in the next diagram, figure 
11, which plots [Oll]AA3727,9/[Om]A5007 against [O m]A5007/H/3. This diagram was one 
of the original diagnostic diagrams suggested by Baldwin, Phillips, & Terlevich (1981) and 
has also been used by Binette et al. (1996) to explore the Am/i model. In our figure both 
the dusty and dust-free models have been plotted, displaying curves of different power-law 
index. In both cases a flatter power-law corresponds to an increase in the [Olll]/H/3 index. 
The observations of Koski (1978) and Allen (1998) are displayed using the same symbols as 
the previous diagram. In this diagnostic, the [Olll]/H/3 provides a standard diagnostic ratio 
which helps to separate the two models, as well as distinguish the different a curves. The 
applicability of the [O Il]/[0 III] ratio as a model diagnostic is clear in this figure, with the 
two model sets separating above logUo = —2.0. The dusty model curves characteristically 
stagnate in ionization parameter above this value, becoming degenerate in Uq in the region 
occupied by the observations. The dust-free models continue past this region to proceed to 
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smaller values of the [0n]/[0lll] ratio. 

In the next diagram, figure 12, we use the [0 n]/[0 ill] ratio to introduce another model 
diagnostic ratio, He IIA4686/H/3. This diagnostic diagram was first introduced by Binette 
et al. (1996) to distinguish the A M /i model, as it was the failings of the standard photoion- 
ization model to reproduce the observed NLR Hell/H/3 ratios that led to the development 
of this model. The Hell/H/3 ratio is distinct in the diagnostic ratios examined so far in 
that it is the dust-free models which stagnate at high ionization parameter, not the dusty 
models. In figure 12 we display curves of different a for both of these models, with a in- 
creasing from right to left. The stagnation of dust-free model in the Hell/H/3 ratio at high 
ionization parameter can be seen in the way the two model sets diverge. The way in which 
the dusty model stagnates in the [O ll]/[0 ill] ratio at high U , while still heading towards 
larger values of the Hell/H/3 ratio, reproduces the observations extremely well, especially in 
comparison to the inverse behavior seen in the dust-free models. The success of the dusty 
model in attaining the high Hell/H/3 ratio with the observed [O n]/[0 ill] values comes from 
its ability to maintain both a high and low-ionization zone at high ionization parameter. 
This is similar to the main idea of the Am/i model, except it is dust which provides the 
absorbed spectrum not "matter bound clouds". The absorption by dust further assists by 
preferentially removing the H ionizing photons to the He ionizing ones, and hence increasing 
the Hell/H/3 ratio. 

Figure 12 is also able to distinguish a, with a large spread in the curves across the 
Hell/H/3 ratio, especially in the dust-free models. Density variations cause similar effects 
to that seen in figure 10, with the [O ll]/[0 ill] ratio providing most of the sensitivity. The 
influence of metallicity is interesting, as it appears to be negligible in the dust-free models, 
yet the dominant parameter in the dusty models at high ionization parameter. 

The consequences of metallicity variations upon the He ll/H/3 ratio are better seen in fig- 
ure 13, which consists of Hell A4686/H/3 versus the temperature sensitive ratio [O lll]A4363/[0 ill] A5007 
(Rom)- The Koski (1978) and Allen (1998) observations are included for comparison. On 
the diagram we display the effects of metallicity variation upon the dusty models. The ef- 
fect of metallicity variation is predominantly seen in the Rom ratio. An increase in the 
metallicity results in a decrease in the nebula temperature and thus a decrease in Rom, as 
metals dominate the cooling processes in nebulae. These affects are discussed in detail in 
Paper 1. In the Hell/H/3 ratio the effects of metallicity are not obvious until Uq > 10~ 2 . 
Although the change in abundance affects the Hell emission, the change in emission with 
metallicity at high U is caused by the change in temperature and the increase in dust. As 
shown in §5.1 of Paper 1, an increase in metallicity leads to drop in temperature in the Hell 
zone relative to the H/3 emitting region, and hence an increase in the recombination to He II 
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relative to H. However this effect only contributes ~10% to the spread. The major effect is 
due to the increase in the dust to gas ratio because of the increase in metal abundance. This 
increases the total dust opacity, which preferentially removes the H ionizing photons to the 
He ones (see figure 2 in Paper 1), thus increases the Hell emitting column relative to the 
Hn column and hence the Hell/H/3 ratio. The spread in model curves is not seen at low 
ionization parameter because the dust is no longer the dominant opacity and hence does not 
control the emitting column. 

The [O m]A4363/[0 ill] A5007 (Roin) ratio is actually one of the standard line ratios, 
used because of its strength and temperature sensitivity. Figure 14 shows a standard di- 
agnostic diagram plotting Roin against [O lll]A5007/H/3. Shown are the effects of density 
variation on both the dusty and dust-free models, with density increasing from left to right 
at log£/o = —3.0. Included with the usual Koski (1978) (crosses) and Allen (1998) (aster- 
isks) observations, are data from Tadhunter, Robinson, & Morganti (1989) (triangles), who 
specifically looked at the "Temperature problem"; the inability of photoionization models 
to reproduce the high Roin ratio. This failure is easily seen in the way the dust-free model 
curves wrap around at high Uo, not reaching the high Rom region occupied by the data. 
What can also be seen in the figure is the success of the dusty models in attaining Roin 
values similar to those observed. This success is due to several reasons, the main being the 
hardening of the spectrum by dust absorption and the contribution of photoelectric heating 
by dust to the temperature. 

Combining two of the previous ratios, the Hell/H/3 versus Roin diagram is one of the 
more powerful diagnostic diagrams. As shown in figure 10 in Evans et al. (1999), it clearly 
separates the different possible excitation mechanisms. It is also a diagram in which the 
standard photoionization model conspicuously fails to reproduce the observations, being too 
cool to get strong Roin and unable to attain the strong Hell/H/3. Figure 15 demonstrates 
this failure, along with the success of the dusty, radiation pressure dominated models in 
reproducing strong Hell/H/5. Curves of different density have been plotted for both models 
sets, with density increasing with increasing Roin- There is a clear separation between the 
dusty and dust-free models, with the dusty models turning over in both ratios before the 
observations are reached. In this instance, the dusty models have too low a Roin ratio 
to successfully reproduce the observations. A flatter power-law and lower metallicity can 
better reproduce the observations due to the increase in temperature and the increase in 
Hell ionizing photons, but such changes disagree with the previous diagnostic diagrams. 
The data itself may also be in error due to the possibility of overestimating such a weak line 
like [Om]A4363 (especially in comparison with a strong one like A5007). Either way this 
diagram indicates that the models still require some further examination. The dust alleviates 
part of the temperature problem, but some additional heating source is still needed, such as 
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small shocks or turbulent heating. 

The Hell/H/3 diagnostic series continues with Hell/H/3 versus [O m]A5007/H/3 in figure 

16. This use of the standard diagnostic ratio provides a clear distinction between the two 
models, plotted here with density the varying parameter. In both models the higher density 
corresponds to a stronger [Olll]/H/3 ratio. This diagram is a simple and clear justification 
of the dusty model over the dust-free one. 

The final diagnostic diagram in the Hell/H/3 series is against [Nn]A6583/Ho; in figure 

17. In 17a we show dusty model curves of different metallicity, along with the Koski (1978) 
and Allen (1998) data. The use of [Nil] /Ha as a metallicity diagnostic has been discussed 
previously but here the Hell/H/3 provides a clear ionization parameter diagnostic. This 
enables a better determination of the metallicity, as well as U . 

As a model diagnostic this diagram is also very good, as demonstrated by figure 17b. 
The variation due to density in both models is negligible when compared to the separation 
of the dusty and dust-free model curves at high ionization parameter (Uq > 10~ 2 ). Even the 
dispersion due to metallicity variations provides only small confusion in distinguishing the 
two photoionization models. At low ionization parameter (Uq < 10~ 2 ), the models become 
indistinguishable as the effects of dust and radiation pressure on dust are relatively small 
below this value. The data from Koski (1978) and Allen (1998) again agree with the dusty 
model. 

The next optical diagnostic diagram examines the partially ionized zone of the NLR 
clouds described in Paper 1 by analyzing the low ionization states of oxygen and nitro- 
gen. Figure 18 plots [N l]A5200/[O ll]AA7318, 24 versus [O l]A6300/[N ll]A6583 for both sets 
of models, with density decreasing from left to right for both model sets. As mentioned 
before, the low ionization species like O I and N I are somewhat untrustworthy in the models 
due to their dependence upon the model termination point, which is at the point where 
H is less than 1% ionized in our models. However the twisting nature of these diagnostic 
curves make this diagram worth investigating. As expected, for logf/o < —2.0 the models 
are very similar, but above this value they diverge. The dust-free model curves turn back 
upon themselves, heading towards small values of the [Nl]/[0n] and [Ol]/[Nll] ratios. The 
dusty models continue to increase in both ratios, largely because the partially ionized region 
continues to grow in the dusty models with Uo, as the dust opacity acts to harden the ionizing 
spectrum. 

The next diagram (figure 19) definitively demonstrates the effectiveness of the dusty 
models. The ratios [S ll]AA6717, 30/Ha versus [N ll]A6583/Ha are plotted for the dusty 
models, displaying the curves of different metallicity. These two ratios together effectively 
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display what was seen with the Veilleux & Osterbrock (1987) curves. When the ionization 
parameter ranges over —4.0 < logC/o < 0.0 the dust-free models cover a range of 3 dex in 
both ratios (as can be seen in figure 1). For the dusty, radiation pressure dominated models 
however, over this large range in ionization parameter each metallicity curve only covers a 
range of ~0.5 dex in both ratios. Even more indicative is that this region of stagnation is 
the same region as occupied by the observations. 

The final optical diagnostic diagram considers the combined ratio [O ll]AA3727, 9x [S ll]AA6717, 30 
/[S ll]AA4067, 76x [O ll]AA7318, 24. In figure 20 the variations in density for the dusty models 
are plotted on this ratio versus the ionization parameter diagnostic [O lll]A5007/H/3. The 
[Oii][Sii]/[Sii][Oii] ratio minimizes the effects of variations in both metallicity and power- 
law index, as the sensitivity of each ratio is cancelled out by its inverse. The ratio has 
both sulfur and oxygen, and short and long wavelength emission lines on both numerator 
and denominator, leaving a ratio which experiences little influence from abundance varia- 
tions and reddening. What it does leaves is a ratio strongly dependent upon density, clearly 
demonstrated in figure 20. 

Though the previous diagrams are only a small sample of the possible optical diagnostic 
diagrams, they provide enough diagnostics such that they can distinguish the excitation 
mechanism and verify that the dusty, radiation pressure dominated model is the correct 
paradigm for photoionization. Together they can also provide estimates for the parameters 
that define the nebulae; density, metallicity and ionizing spectral energy distribution. 

3.3.3. Near-Infrared Diagnostic Diagrams 

Though the Near- Infrared experience little extinction due to dust, there is a scarcity 
of strong lines in this spectral region and thus limited choices for diagnostics. The noble 
elements, Ne and Ar dominate this spectral region and provide good diagnostic emission 
lines as they are not depleted onto dust. 

The first diagnostic, figure 21, 
plots [Ne Vl]A7.652/zm/[Nell]A12.8/im versus [Ar vi]A4.53/zm/[Ar n]A6.98/im and displays the 
effects of variation of a upon the dusty model. Note that as they both contain emission lines 
from high ionization species, the ratios from the model rapidly disappear as the ionization 
parameter becomes small (Uo < 10~ 2 ). The IR diagnostic diagram is similar to both the 
UV and optical diagrams in that it demonstrates the stagnation of the dusty models at 
high ionization parameter. It also shows the same divergence of the dusty and dust-free 
models above \ogU = —2.0, though not indicated here. This particular diagnostic diagram 
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is interesting in that it displays a strong sensitivity to variations in a, with both the influence 
of metallicity and density being negligible. The sensitivity to metallicity is removed through 
having ratios of the same element, making this diagram a very good diagnostic of the ionizing 
SED. The sensitivity to a arises from the difference in ionization states of the numerator 
and denominator in both ratios. 

The next diagram also displays a similar sole sensitivity to a. 
Figure 22 plots [NeVi]A7.652//m/[Neii]A12.8/mi versus [Ar iii]A8.98/im/[Ar ll]A6.98//m dis- 
playing curves of different power-law index for both the dusty and dust-free models. Ap- 
parent in the diagram is the divergence of the two model sets at large values of ionization 
parameter (C/ > 10~ 2 ). Also conspicuous is the characteristic stagnation in the dusty models 
at these high values of ionization parameter. The power-law index of each curve decreases 
with increasing value of the Argon ratio for both sets of models, with —2.0 at the top to —1.2 
for the bottom curve. The use of this diagram as an a diagnostic is clear, and the negligible 
reddening in this spectral region make this an ideal diagnostic for heavily obscured NLRs. 
Note that the Ne ratio has a large separation at low Uq, which is not obvious due to the 
different scales on each axis. This is as when the power law steepens, it drops the number 
of Ne V ionizing photons below a critical point and Ne VI becomes very weak. 

With these two diagrams we have what appears to be rare in the UV and optical; a 
pure power-law diagnostic. Thus, if an active galaxies emission lines are measured for all 
three spectral regions, the use of line diagnostic diagrams can not only put strong constraints 
on the excitation mechanism but also the density, metallicity and ionizing spectrum of the 
emission line region. 

4. Discussion 

Through the diagnostic diagrams and comparisons with observations in the previous 
section we have demonstrated the validity of the dusty, radiation pressure dominated models 
as the paradigm for photoionization. The ability to reproduce both the observations and 
observed clustering of narrow line regions with similar parameters in diagnostic diagrams is 
strong evidence for the model. However, this work does not negate all work based upon the 
standard photoionization model. The active galaxies which were believed to be photoionized 
still appear to be so. The parameters estimated for these regions by the standard model 
are in general still correct. What this new paradigm does provide is a better diagnostic for 
the excitation mechanism, an explanation of how and where the standard model failed and 
better diagnostics to obtain more accurate estimates for the parameters that define these 
regions. 



-23- 



Neither does this new paradigm prove wrong the previous improvements upon the stan- 
dard model, such as the Am/i model (Binette et al. 1996, 1997) and the "locally optimal 
emitting cloud" model (Baldwin, Ferland, Korista, & Verner 1995; Ferguson et al. 1997). 
Rather, this work, itself an improvement upon the standard model, provides a self consis- 
tent and self-contained physical basis for the emission of both high- and low-ionization lines 
without the necessary combining of clouds of different densities and ionization states. 

As such, the dusty, radiation pressure dominated models can be considered as an expla- 
nation for the assumptions of the previous models. For example, the structure of the clouds 
in the dusty model, a depiction of which is displayed in figure 1 of Dopita et al. (2002), closely 
matches one of the geometrical distributions suggested for the A M /i models (Binette et al. 
1996, figure 4b). The dusty model provides both the low ionization, ionization bounded com- 
ponent and the high ionization, matter bounded component within a single model, as shown 
in the ionization structure diagrams given in Paper 1. It does this through self-shielding, 
such that the high-ionization component sees the direct spectrum while the low ionization 
component near the ionization front sees the self absorbed spectrum. The dusty models 
also provide an explanation for the pressure difference between the matter bounded and 
ionization bounded components of the Am/i model through the effects of radiation pressure 
on dust. The final supporting point is that the dusty, radiation pressure dominated models 
also answers the problem emission line ratios that the A M /i models was originally created to 
solve. As shown in figure 13 and 15, the dusty model is able to produce strong Hell/H/9 and 
strong [O lll]A4363/[0 ill] A5007 (Rom), through the combination of ionization zones and 
dust photoelectric heating. The main distinction between the results of the dusty models 
and the Am/i models is also the major improvement; the stagnation in ionization parameter 
characteristic to the dusty model. 

Just as the dusty, radiation pressure model does not prove incorrect the previous pho- 
toionization models, neither does it prove wrong the shock excitation models. The shock 
models and the photoionizing shock+precursor models require some form of mechanical en- 
ergy input. In the absence of such an input the emission line clouds must be photoionized. 
However, if a jet or strong outflow is present within the AGN then shock excitation will play 
some part. In most cases the contribution from shock excited clouds is either indistinguish- 
able or negligible compared the contribution from photoionized clouds. However in some 
cases, the emission from the shock excited clouds can become noticeable or even dominate 
the emission from photoionized clouds. Such a case is seen in the UV diagnostic shown in 
figure 5 when compared to the shock models and observations from figure 14 in Inskip et al. 
(2002a). Most of the radio galaxies are reproduced quite well by the dusty models. However 
a few of the observed high-z radio galaxies have a lower Clll]/Cll] ratio then is obtainable 
for the values of [Nell]/[Nev] seen. However, as discussed in Inskip et al. (2002a), these 
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outlying galaxies can either be easily fit by the shock+precursor model or a combination of 
shock+precursor and dusty model. Inskip et al. (2002a) also suggest that the balance be- 
tween the two components depends upon the radio source size and hence jet structure, which 
may indicate some form of jet-cloud interaction. Even so, these observations demonstrate 
that more needs to be understood about how these two methods of excitation can coexist 
and how they interact. 

5. Conclusion 

First introduced in Dopita et al. (2002) the dusty, radiation pressure dominated pho- 
toionization model provides a self-consistent explanation for the emission from the narrow 
line regions of AGN. Within this work we have continued the examination of this new 
paradigm begun in Groves, Dopita & Sutherland (2003). 

Through the comparison of observations on standard optical line diagnostic diagrams 
of Veilleux & Osterbrock (1987) and the UV diagnostics of Allen, Dopita, & Tsvetanov 
(1998) we have clearly demonstrated the validity of the dusty model as the new paradigm 
for photoionization. The dusty model, through the stagnation of the ionization parameter 
at large values, provides a simple explanation for the small variation of observed Seyfert 
NLR ratios. This stagnation is due to the effects of radiation pressure upon dust and is 
characteristic to these models. The significant point is that the dusty model is able to do 
this over both optical and UV ratios, without depending upon large variations in other 
parameters such as density or metallicity. 

As well as verifying the dusty model,we have also explored the effects that the variation 
of density, metallicity and ionizing spectrum have upon both the new dusty paradigm and the 
standard photoionization models. With this exploration, we have demonstrated that several 
known line ratio diagrams can be possibly used as diagnostics of these three parameters. 

In addition to the previously explored diagnostics, we have introduced several new line 
diagnostic diagrams, covering UV, optical and IR ratios. These diagnostics provide further 
tests for the dusty model as well as providing diagnostics for metallicity, density, ionizing 
spectrum and ionization parameter. 

These results not only provide an explanation for what has not been a fully understood 
observation for years but also provide ways in which to understand further the processes 
involved in the NLR and extended NLR of AGN. 
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Table 1: Solar Abundance Set 



Element 


Abundance 21 


Element 


Abundance* 


H 


0.000 


Al 


-5.51 


He 


-0.987 


Si 


-4.49 


C 


-3.61 


s 


-4.80 


N 


-4.20 


CI 


-6.72 





-3.31 


Ar 


-5.60 


Ne 


-3.92 


Ca 


-5.65 


Na 


-5.68 


Fe 


-4.54 


Mg 


-4.42 


Ni 


-5.75 



a All abundances are logarithmic with respect to Hydrogen 
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Table 2. Depletion Factors 



Element Depletion' 



H 


0.00 


He 


0.00 


C 


-0.30 


N 


-0.22 





-0.22 


Ne 


0.00 


Na 


-0.60 


Mg 


-0.70 


Al 


-1.70 


Si 


-1.00 


s 


0.00 


CI 


-0.30 


Ar 


0.00 


Ca 


-2.52 


Fe 


-2.00 


Ni 


-1.40 



Depletion given as log(X/H) gas - log(X/H) ISM 
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Fig. 1. — The standard optical line diagnostic diagrams of Veilleux & Osterbrock 
(1987) (V&O) showing [Nll]A6583/Ha, [Ol]A6300/Ha and [S ll]AA6717, 30/Ha versus 
[Olll]A5007/H/3. The data points are from (V087) and curves are as labelled in the keys. 
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Fig. 2.— UV diagnostics from Villar-Martin, Tadhunter, & Clark (1997, VTC97) and Allen, 
Dopita, & Tsvetanov (1998, ADT). Observations are from ADT, with asterisks represent 
high-z radio galaxies and triangles representing HST observations of three Seyfert 2 NLRs. 
Dust-free curves are marked with diamonds and dusty curves are marked with squares. Arrow 
represents magnitude and direction of 3 visual magnitudes extinction of models through a 
dusty screen. 
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Fig. 3. — UV diagnostic from ADT, demonstrating the effects of metallicity on the dusty 
model (a) and differences between the two models (b). a decreases from top to bottom for 
both dusty (square) and dust-free (diamond) model curves in (b). Observations (triangles) 
are from ADT. Reddening arrow as in figure 2a 

Fig. 4. — UV diagnostic from ADT. Triangles represent HST observations from ADT, aster- 
isks represent observations of Seyfert 2s from Koski (f 978) and crosses represent observations 
of Seyfert 2s from Allen (1998). In (b) [Olll]/H/3 increases with a for both model sets. 
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Fig. 5. — UV diagnostic diagram from Best, Rottgering, & Longair (2000b); Inskip et al. 
(2002a). Observations from Inskip et al. (2002a) with asterisks representing 3C galaxies and 
triangles 6C galaxies. In (b) a decreases with increasing Ciii]/Cii]. Dusty model curves are 
marked with squares and dust-free models with diamonds. 
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Fig. 6. — UV diagnostic diagram of a dusty model with Z = 2Z Q and a — —1.4 as marked, 
showing curves of different density as labelled. Lines of constant ionization parameter are 
as labelled. 

Fig. 7. — UV diagnostic diagram. Labelled as in previous diagram. This diagram demon- 
strates both a sensitivity to metallicity, due to the component of secondary Nitrogen, and 
a sensitivity to ionization parameter, revealing the stagnation due to radiation pressure on 
dust. 

Fig. 8. — UV diagnostic diagram. Labelled as in previous diagram. Shows a strong metal- 
licity sensitivity on both axis due to Nitrogen in each ration. 

Fig. 9. — UV diagnostic diagram. Labelled as in previous diagram. Both metallicity and U 
sensitive as in figure 7. 
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Fig. 10. — Near-UV optical diagnostic diagram showing density variations on both the 
dusty (squares) and dust-free (diamonds) models. Lines of constant are Uo marked on both. 
Density increases from top to bottom. Observations from Koski (1978) (asterisks) and Allen 
(1998) (crosses). Reddening arrow indicates magnitude and direction of 3 Ay extinction on 
models. 

Fig. 11. — Optical diagnostic diagram suggested by Baldwin, Phillips, & Terlevich (1981) 
showing a variation in dusty and dust-free model. Curves and data marked as in previous 
diagram. Reddening arrow indicates 10Ay extinction here, a increases with [Olll]/H/3. 

Fig. 12. — Binette et al. (1996) diagnostic diagram showing a variation in dusty and dust- 
free models. Models and data marked as in previous diagram, a increases from right to 
left. 

Fig. 13. — Diagnostic diagram showing metallicity variations of dusty model. Curves and 
data marked as before. 
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Fig. 14. — Temperature sensitive diagnostic diagram showing density variation of both dusty 
and dust-free model curves. Density increases left to right at U = 1CT 3 . In addition to the 
Koski (1978) (asterisks) and Allen (1998) (crosses) data, data from Tadhunter, Robinson, & 
Morganti (1989) is shown. 

Fig. 15. — Helium - Oxygen diagnostic diagram showing density variations on dusty and 
dust free models. Density decreases from top to bottom for both sets of model curves. 

Fig. 16. — Helium - oxygen diagnostic diagram with density variations of both models. 
Density increasing with increasing [Olll]/H/5. Symbols as before. 
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Fig. 17. — Helium - nitrogen diagnostic diagram showing metallicity variation of dusty model 
(a) and density variations of dusty and dust-free models (b). Symbols as before. 
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Fig. 18. — Density variations of both models in low ionization species diagnostic diagram. 
Density increases right to left for both sets of curves. 

Fig. 19. — Metallicity variation of dusty model on nitrogen - sulfur diagnostic diagram. 

Fig. 20.- [On]AA3727,9x[Sn]AA6717,30 /[Sii]AA4067,76x[Oii]AA7318,24 versus 
[O lll]A5007/H/3 diagnostic diagram showing density variation of dusty model. 
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Fig. 21. — Infrared diagnostic diagram showing curves of different a for the dusty model. 

Fig. 22. — Comparison of the dusty (square) and dust-free (diamond) models in the Infrared 
diagnostic diagram. Curves of different a are shown, with a increasing from top to bottom 
for both models. Note the different scales on each axis. 



